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ABSTRACT 


An analytical investigation has been made into the stability 
of the burninq configuration of a single coaxial injector surrounded 
by similar injectors. The stability criteria is based on an average 
pressure difference along the boundaries of the adjacent stream tubes 
as calculated using Spaulding's numerical method. 

The results indicate qualitatively that there is a tendency for 
the injectors to have different burning configurations. It is believed 
that the configuration achieved is random; however once the burning 
configuration is established, it is believed to persist. These results 
are consistent with previous experimental observations. 
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I. INTRODUCTION 


Over the nast several years a series of experiments have been 
(1 2 3 4 ) 

conducted tor the purpose ot determining the turbulent 

mixing parameters in a liquid rocket engine. A tracer gas diffusion 
method with helium gas injected at a point along the combustion chamber 
centerline was used. (See figure 1) Downstream of the helium injector, 
gas samples were withdrawn at a series of points along a major diameter. 
Later these samples were analyzed to determine the helium concentration 
in each sample. The turbulent mixing parameters were calculated using 
the theories of turbulent diffusion and the measured helium concentration 
values. 

The same small rocket engine was used for all the tests. This en- 
gine had a cylindrical combustion chamber S.94 cm in diameter and approx- 
imately 37 cm long with a 2.S4 cm throat. The initial series of runs 
we**e made using liguid oxygen and heptane as propellants and a like-on- 
like imoengement injector. During these runs the sampling and measuring 
techniques were developed, and consistent and repeatable data were ob- 
tained. 

A series of runs were then conducted using gaseous oxygen and gaseous 
hydrogen and later liguid oxygen and gaseous hydrogen as propellants, 
for both of these propellant combinations, a seven element coaxial in- 
jector was used as shown in figure 1. 

Throughout both series of hydrogen-oxygen runs, a large amount of 
scatter was observed in the helium concentration measurement both within 
a run and between runs. Six gas samples could be taken per run, and it 
became customary to make four runs at each test configuration thus having 
twenty-four data points to define a helium concentration curve. 
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Propellant Sample 
Injector Lines 
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Mass Flow Rate GOX = 0 ,?] «g/sec 
Mass Flow Rate CH^ = p,0\R Rg/sfc 
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FIGURF 1 - EXPER I MFNTAL CONFIGURATION 



The six ports in the sample probe were distributed along the 
entire diameterand the probe could be moved along the diameter in 
small fixed increments between runs. Thus subsequent runs would 
provide data points to fill in the spaces between points from the 
previous runs. Figure 2 shows the data from a typical set of four 
runs. The large amount of scatter is readily apoarent. 

Repeated efforts were made to determine the cause of these in- 
consistencies, All measuring and metering equipment was recalibrated 
and/or replaced. The entire system was inspected for leaks many times. 
An error analysis which indicated the accuracy of the helium concentra- 
tion measurements should have been within about ± 5% was made; however, 
the actual scatter was many times this value. 

It was finally concluded that Lne concentration measurements must 
be representative of the conditions in the engine. The possibility 
that the problem was caused by some tyDe of combustion instability in 
the usual sense was ruled out. Even low-frequency combustion insta- 
bility will be in the order of several hundred hertz. The sampling 
time for each of these runs was approximately 1.5 seconds. Certainly 
variations caused by combustion instability would be averaged in 1.5 
seconds. If it were some type of very low frequence "chugging," the 
existance of the condition should be apparent on either the chamber 
pressure or thrust measurements; however, it was not. 

It is the purpose of the work herein to investigate analytically 
the possibility that the engine had more than one stable mode of 
burning. 
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To analyze such a complex flow situation involving a bi-orooellent 

system with turbulent combustion is very difficult even under steady 

state conditions, and to include the time variable would be out of the 

question. The approach chosen was the Spaulding steady state computer 

( 5 ) 

analysis given by Gosman et al. This method was used to calculate 
the pressure along the boundary of a single element stream tube under 
different conditions of stream tube diameter and propellant flow rate. 
The results of these calculations would then be analyzed. 

It was also decided that due to the nature of the problem and the 
large amount of computer time required for each case considered a 
very simplistic approach would be taken. Only gross trends would be 
sought rather than numerical bounds. 


CALCULATION Of PRESSURE. COEFFICIENTS 

A. The Spaulding NumericaJ Method 

The experimental combustion chamber had seven coaxial propellant 
injectors, thus the flow was three-dimensional . Mathematical lv this 
is a very difficult problem. To simplify the problem, it was assumed 
that each injector formed an axisymmetric jet with turbulent mixinq 
and chemical reaction. The axi sv mmetric jet was assumed to be confined 
in a constant diameter cylindrical stream tube. This boundary condition 
is believed to closely approximate the case of a simile injector surrounde 
bv similar injectors. 

The Spauldimi numerical method is applicable to this situation and 
was used to calculate the pressure dist ibution alone the boundary. The 
uniqueness of this method, which is restricted to two-dimensional or 
axisymmetric problems, lies in the transformat ion of each of the conser- 
vation equations into a standard form which is a nonlinear partial differ- 
ential equation of the elliptic type. The standard equation is solved bv 
a finite difference procedure tisine the iiuuss-Seidel iterative technique. 

A very important consideration ; s the method of model inq the turbulen 

flow. In the present calculat ions . the turbulent flow was modeled bv an 

equation of the effective viscosity, ’’he equation was that proposed bv 
. (b) 

bos man et al. 

The details or this calculation procedure are qiven by desman et al 
includinq a complete listinq of the proqram. Tou presents details 
of the calculat i nq procedure as applicable to the conf iqurat ion considere 


herein . 


B . Con figu ration considered 

The configuration considered was a single injector element with 
a circular stream tube of constant radius. The same injector dimensions 
we- e maintained throughout and are those shown by figure 1. The nominal 
or reference case was assumed to be that in which the cross sectional 
area of the stream tube is one seventh of the chamber area. Under this 
assumption, the stream tube radius is 1.123 cm. a ree longitudinal 
boundary was assumed, i.e., the gradient of the longitudinal velocity 
with respect to the radius equals zero. The nominal or reference in- 
jector velocities for oxygen and hydrogen were 233.5 m/sec. and 780 m/sec. 
respectively. 

Computer runs were made for f ve cases as listed in TABLE I. 

These were: nominal, radius increased 10"., radius decreased 10 , radius 

increased and injector velocities decreased by 10 , and radius decreased 
and injector velocities increased by 10". 

The longitudinal boundary and injector face pressure coefficients 
as calculated by the Spaulding computer program are plotted on figure 3 
and 4, respectively. The pressure coefficient used is referenced to the 
center line pressure, P r , far down stream (z=36.88 cm) and defined as 

C = P/P - 1 (1) 

p r 
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FIGURE 3 PRESSURE COEFFICIENT ALONG LONGITUDINAL BOUNDARY 
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TABLE I -Listing of Computer Cases 


Case ,\R/R 0 aV/V 0 *(») V ox^ n,/stn '' V fu^"’ s0 


* 

1 

0 

0 

1.123 

233.5 

780.0 

2 

+ 0.1 

0 

1.235 

233.5 

780.0 

3 

-0.1 

0 

1.011 

2 33.5 

780.0 

4 

+ 0.1 

-0.1 

1 .235 

210. 1 

702.0 

5 

-0.1 

+ 0.1 

1.011 

2S6.S 

858.0 
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III. Kl SUITS 


A. Pevelopmont of Pressure function 

A s mont ioned in t ho iotroduct ion. it is believed th.it the complovi tv 
of this problem ami the limitations of tho analytical tools at ha ini make 
only the simplist analysis .iust i t iable. 1 Inis . in tho analysis which follows, 
ijross approsimat ions ami simplifying assumptions aro made with tho umior 
stamiimi that possible trends only art' to bo implied from tin' results 
obt a i nod . 

Inspection of tin' curves of injure ' shows th.it cliamn's in stream 
tube radius ami injector velocities have an effect on t he lono i t tniina 1 
boundary pressure distribution in the rovjion ’ ran about ont' to si\ centi 
meters downstream of the injector, further downstream, t ho effect is very 
small. Ihis is believed significant because i* tin' stream tube radius is 
altered by a pressure difference in tin* upstream region, the ihatnje will 
persist in the downstream tea ion. This would »jive some just i t icat ion tor 


assumino a i onstant diameter stream tube. 

fioure -1 shows that chamjes in stri'am tube radius and injector ve’.n : c. 
have a latnje effect on t he injector face pressure ioe''uient ’hose itve* !cr 
favi' pressure variations are laroo enouuh to cause significant chamjes in the 
pressure across the injector and thus influence the propellant 'low rates 
The pressure data hue hern referenced to tin' pressure tar downstream 
because it is assumed that tins pressure will be tin' controllimj boundary 


condition when several injectors are adjacent tv' one another. 

1 tie chamjos from nominal in tin' pressure distribution caused bv stream 

tube radius changes and injector velocity chamjos are also shown on future : 

and 4. Average chamjes in the lomji tudina 1 boundary pressure coefficient. 

\f . and average chamjes in inioctor face pressure coefficient. U\, . have 
h f 

t'een assumed as indicated by tin' straight lines on the oraohs . 
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Those average changes were approximated bv oguations as follows: 

ac p o.12|ar/r 0 | - 0.2 (aV'v 0 ) (r) 

b 

and 

Af ( , 0.7 (aR/R 0 ) ♦ 0.4 (aV'V 0 ) (3) 

The accuracies of these oguations are illustrated by figure b. 

1 he relationship between the pressure drop across the injector and 
the injection velocity can be approximated bv the incompressible flow 
eguat ion. 

,\P. • (fl P ♦ 1 ) eV‘ /: (i ' 

i nj 

where \P . . \P . . ♦ \(\P. .). Assuming that variations in t 1 ow rate 

' inj inj ^ i n.i 

do not significantly change tin' upstream injector manifold pressures, then. 

\l\P. .) -P (ac ). Pi fferent iat ing eguat ion 4 for small deviations 
' ' in.i r P • 

from nominal conditions it follows that: 

At',, - 1 T. (fl P - 1) eV 0 *' (A V V 0 ) ^ 

Substituting numerical values tor the t'xper iment a 1 cent igurat ion . the 
relationship can be represented to within IP tor both oxvgen and hvdrogon 
streams bv the eguat ion 

U'p -1.1 (AV V 0 ) (b) 

Combining oguations j and p gives an expression tor the average boundarv 
pressure coefficient in terms of the change in radius. 

Af p 0.1.' ’ \R R A ♦ P. IP ( \R R ,1 C 1 

b 

Consider now a cent igurat ion where one central injector element is 
surrounded by N other injectors and when the stream tube of the central 
injector decreases in radius the stream tubes ot the N surrounding in- 
jectors increase in radii, and visa versa. Assume that the cross section 


ot all N ♦ 1 stream tubes remains constant as it would in a small engine. 
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FIGURE 5 ACCURACY OF EMPERICAL EQUATIONS FOR 
VARIATION IN PRESSURE COEFFICIENT 


Define a pressure parameter 


* = AC p 


- AC, 


( 8 ) 


b ) o h ) s 

to represent the pressure difference between tin' central and surrounding 
stream tubes. The subscripts , and s refer to larqe and small, respectively 
A positive value of | indicates an ustable condition because tin' larqe 
stream tube will tend to become larqer. 

Under the conditions outlined above, if the central stream tube is smalle 
than the surroundinq stream tubes, i.e.. \R, 0 

* - -n. l ?|.\R C /R I (l - I NI - n.io (aR c /R o ) (l ♦ l/Nl (d) 

and if the central stream tube is larqer, i.e., \R s 0 

1 . 1 2 1 AR ,/Rj (1 - 1/N) ♦ 0.10 (AR /R o ) (1 ♦ 1/N ) (18) 


where \R . is chanqe in radius ot the central stream tube. 
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B . Di s cussi o n of Results 

The pressure function, is shown on figure 6 as a function 
of the fractional change in the radius, aR c /R q , and the number of surround- 
ing injectors, N. It must remember that the pressure function represents 
the pressure coefficient of the larger stream tube minus the pressure 
coefficient of the smaller no matter which happens to be in the center. 
Therefore, a positive value of pressure function, 4 -, means that a larger 
stream tube will tend to get larger, i.e., there is a destabilizing effect. 

The pressure for a two iniector combination is represented by N = 1 
and is seen to have a very large destabilizing tendency. According to this 
analysis, the two injector combination would always tend to burn in an 
asymmetric configuration. 

The case of a single injector surrounded by six other injectors where 

the six surrounding injectors have identical burning patterns is represented 

by N = 6 . For aR /R --0, i.e., the center stream tube is smaller, the system 
co 

is seen to be neutrally stable. On the other hand for aR c /R q > 0, i.e., the 
center stream tube is larger, the results indicate a large destabilizing 
influence. 

While this analysis is admittedly very crude, it does give an indication 
of the type burning patterns which might tend to form. In particular, it 
indicates that the case in which all injectors have the same burning pattern 
is not very likely. Figure 7 shows some of the burning patterns which it is 
believed will tend to form. The particular pattern which actually forms mav 
be a matter of chance depending on which injector ignites first. 

Likewise, the crudeness of the analysis makes it impossible to make 
any predictions regarding the details of the final steady state burning con- 
figurations. Certainly, the steady state configuration would be decidelv 
three-dimensional and many non-linearities would be present. 
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The destabilizing influence Dredicted above is for a particular 
coupling between the propellant manifold and the combustion chamber and 
for a constant manifold pressure. Both of these situations are believed 
to be reoresentati ve of the typical hydrogren oxygen engine; and while 
the numerical values of the pressure function might change, the trends 
would be the same. 

These results do not prove that the scatter in the experimental data 
discussed in the introduction was the result of the stable but dissimilar 
burning of the individual injectors. It does, however, indicate that such 
burning might be possible and shows the mechanism whereby it could occur. 

The analysis indicates that the injector face pressures surrounding 
the individual injectors could differ by an easily measurable amount. If 
such measurements were made, they should indicate whether or not the dis- 
similar burning condition actually exists. Also, if such a condition were 
found, it would add significant creditability to the Spaulding method. 
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IV. CONCLUSIONS & RECOMMENDATIONS 


The conclusions drawn from this work are: 

(1) The Spauldinq numerical method can be used to predict that 
the individual injectors of a mul ti -coaxial injector system 
will tend toward dissimilar burning confiqurations. 

(2) There may be several stable steady-state burninq confiqurations, 
and the particular configuration occurring may depend on the 
iqnition process. 

(3) The dissimilar burning condition is a result of couplinq with 
the propellant manifolds. 

(4) The dissimilar burning conditions could be the cause of the 
larqe scatter observed in previous experimental turbulent 
diffusion data. 

It is recommended that: 

(1) Further analysis should be conducted in order to define the 
steady-state burninq configuration(s) of multi-coaxial injector 
systems . 

(2) Experimental verification of the destabilizing tendency predicted 
herein should be souqht roth for the purpose of determining the 
burning characteristics and for tht purpose of evaluating 

the Spaulding numerical method. 


18 


